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Abstract: In this paper, we study an extension of the standard model with a vector-like 
generation of leptons. This model provides a viable dark matter candidate and a possibility 
to enhance the Higgs decay rate into a pair of photons. We evaluate constraints from 
electroweak precision tests and from vacuum stability, and find that the latter provide an 
upper limit on the lepton Yukawa couplings. A large enhancement of the Higgs di-photon 
rate can therefore only be obtained when the mass of the lightest charged lepton is close to 
the LEP limit. The relic density constraint suggests a co- annihilation scenario with a small 
mass difference between the lightest charged and neutral leptons, which also weakens the 
LEP limit on the lightest charged lepton mass and allows for larger Higgs di-photon decay 
rates. Cross sections for direct detection of the dark matter candidate are calculated, and 
prospects for detecting the new particles at the LHC are discussed briefly. 



1 Introduction 



The discovery of a scalar resonance consistent with the SM Higgs boson [1, 2] marks a great 
success in the early run of the CERN Large Hadron Collider (LHC) and its associated 
experiments. It is now conceivable that electroweak symmetry is indeed broken by the 
Higgs mechanism, namely through the vacuum expectation value of a fundamental Higgs 
doublet. With a Higgs mass of about 125 GeV, this model is fully renormalizable and 
potentially stable up to very high scales. 

However some mysteries remain that suggest additional particles being present at the 
electroweak scale. First, the existence of dark matter has been established through astro- 
physical observations beyond reasonable doubt, however its nature and properties remain 
largely unknown. With all known matter content in the SM being related to the weak 
scale, it is at least a reasonable assumption that also the dark matter particle is connected 
to the electroweak scale. 

The second mystery is provided by the Higgs boson itself. While the existence of a new 
resonance has very recently been established, it remains to be seen whether its properties, 
in particular the production cross sections and decay branching ratios, agree with the very 
precise predictions from the SM. The decay channels that are most sensitive to new physics 
effects are the loop induced decays of the Higgs to pairs of photons or to a photon and a 
Z boson. And indeed the rate of di-photon events shows the largest and most consistent 
deviation from SM expectations, with both 2011 [3, 4] and 2012 [5, 6] data from both the 
Atlas and CMS experiments suggesting an enhancement of the branching ratio of about 
50% (but still within SM expectations at the two sigma level 1 ). 

Motivated by the above, here we study a simple extension of the SM that provides 
a dark matter candidate and allows for an enhancement of the Higgs branching ratio to 
photon pairs. The model introduces a vector-like fourth generation of leptons, namely 
SU(2) doublets t = (^l^r)' charged SU(2) singlets e = (e^e^) and neutral singlets 
v = (z/p, In the limit where the vector-like masses vanish the particle content is that 
of a fourth generation of leptons (indicated by a single prime) and an exact copy with 
opposite chirality (double primed). 

The remainder of the paper is organized as follows. In the next section, the particle 
content of the model is introduced and the mixing in the charged and neutral sector is 
discussed. Electroweak precision tests and other constraints on the model are presented in 
Sec. 3 while the modified Higgs boson properties are analyzed in detail in Sec. 4. In Sec. 5 
constraints on the magnitude of the Yukawa couplings from vacuum stability considerations 
are derived. Dark matter properties are explored in Sec. 6 before we conclude in Sec. 7. 
Formulas for the electroweak precision observables are presented in detail in the appendix. 

Very recent analyses of the combined ATLAS and CMS data sets suggest that the deviation of the 
di-photon signal from SM expectations is at the 2<j-2.5cr level [7-11]. 
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u' R 


l' R 


p" 


< 


Quantum Numbers 


(1,2,-1/2) 


(1,1,-1) 


(1,1,0) 


(1,2,-1/2) 


(1,1,-1) 


(1,1,0) 



Table 1. Particle content of the model with SU(3) x SU(2) x U(l) quantum numbers, in the 
chirality basis. 

2 The Model 

The model introduces one additional family of standard model (SM) like leptons consisting 
of a left-handed SU(2) doublet £' L and the corresponding charged and neutral singlets 
and v' R to allow mass generation through the Higgs mechanism. In addition one mirror 
family of leptons is introduced with identical quantum numbers, but opposite chirality. This 
makes the model manifestly anomaly free, since the fields of opposite chirality combine into 
vector-like multiplets. The full particle content introduced in this model is displayed in 
Tab. 1. 

The Lagrangian density of the model is given by 

1 i 

C = -m^t'i - m e e'le' R - m v v'{y' R - -m'u R c u R - -rri'v^v^ 

- Y^(F L H)e' R - YMrH^)u R - Y?(&B)£ ~ Y^rH^v'l + h.c. (2.1) 

where the standard kinetic terms are omitted. We have written down all mass terms and 
couplings that are allowed by gauge invariance, including Majorana masses for the neutral 
singlet neutrinos. Couplings that would mix the new leptons with SM leptons are set to 
zero for the remainder of this paper. Vanishing of these couplings can be guaranteed by 
imposing a parity symmetry on the new sector, similar to [13]. For the remainder of this 
paper, we will always assume that such a symmetry is present. 

In terms of vector- like multiplets, the particle content can be written as £ = (£' L , £ R ) and 
e = (e^,e R ). To avoid confusion between the flavor and mass eigenstates, the notational 
conventions are chosen as follows: 

• Flavor basis: Lower case particle names v, e and mass parameters m x . 

• Mass basis: Upper case particle names .^1,2,3,4 and £1,2, and corresponding masses 
M Ni , M E .. 

There are two interesting limits of the above Lagrangian that can be realized by im- 
posing a discrete symmetry on some of the new fields. First, one can impose a second 
parity symmetry under which the new SU(2) singlet fields are odd, while all other fields 
are even. This symmetry forbids Yukawa couplings between the new fields such that 
Y' c = Y" = Y£ = Y" = 0, and the masses for all fields come exclusively from explicit 
vector-like mass terms in the Lagrangian, and are not affected by electroweak symmetry 
breaking. In this limit the new fields can be decoupled from the SM without any observable 
low energy effects. While we will not discuss this possibility in detail, the decoupling limit 
is useful for checking the calculation of the electroweak S and T parameters. 
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Alternatively, one can impose an additional discrete symmetry under which the mirror 
leptons are odd, such that the explicit mass terms mg, m e and m v are forbidden. In this 
limit, the mirror sector can not mix with SM leptons, however the new SM like leptons will 
behave like an ordinary fourth lepton family, with corresponding experimental signatures 
and limits. A model with a similar leptonic sector was also considered in [14, 15]. For 
the remainder of this paper we will focus on the most general case, where only the parity 
symmetry that forbids mixing with ordinary SM leptons is present. 

In the limit where the explicit mass terms vanish, the spectrum of the model can 
easily be derived from the Lagrangian. After electroweak symmetry breaking, there are 
two charged leptons with masses Y' c v and Y"v, where v = 174 GeV is the Higgs vacuum 
expectation value (VEV). In the neutral sector the two massive neutrino states are further 
split when the Majorana masses are nonzero, such that there are four neutrinos with masses 

% 3 =/f+^4' (2 - 2) 

Mn sa = ^ + Y"^± y . (2.3) 

When all masses and Yukawa couplings are non-zero, the spectrum is slightly more 
complicated, since now there is mixing between the ordinary and the mirror leptons. In 
the charged sector, the 2x2 dimensional mass matrix M c is defined as 



C D ( e' L e'l ) M c 



R 



+ h.c. with M t 



Y' c v mg 
m P Y"v 



In the neutrino sector all four states are mixed with a symmetric mass matrix Ai, 
by 



\ {< v'r < M n 



( Y'v m e 



+ h.c. 



with 



Mr 



' ■ \ 

Y^v m' m„ 
m e Y"v 



V m v Yl' 



v m 



The mass matrices can be diagonalized as follows: 

M cD = U{M C U R , 
M nD = V T M n V, 



(2.4) 



given 



(2.5) 



(2.6) 
(2.7) 



where M c d an d M n D are diagonal positive semi-definite matrices and Ul, Ur and V are 
unitary. We denote the mass eigenstates as E\ and Ei in the charged sector and iVi_4 in 
the neutral sector with corresponding masses Mg. and Mjv i} ordered such that Mj < Mj 
when i < j. 
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Figure 1. Decoupling of effects on electroweak precision observables. The green (grey) shaded 
region shows the allowed parameter space as function of the mass splitting parameters Am' = 
{Y' c — Y^)v and Am" = (Y" — Y")v. Also shown are contours of constant AT (solid) and AS 
(dashed). On the left the vector-like masses are all vanishing, while the right plot is for mi — m e — 
m v = 400 GeV. 

3 Electroweak Precision Tests and Direct Limits 

In the limit where new particles are heavier than Mz/2, the effects on electroweak precision 
observables can be very well estimated by calculating the contributions to the so called S 
and T parameters [16]. 

We have calculated the contributions to the oblique electroweak parameters S and T 
for the most general allowed mass matrices M c and M n . The formulas are given in the 
appendix. To our knowledge, these results are not available in the literature. Formulas 
for limiting cases are however available, and have been used to verify our results. In 
particular, the case of a chiral generation of leptons that mixes with the SM was treated 
in detail in [17], while some results for mixed vector-like fermions can be found in [18]. 

For a Higgs mass of mn = 125 GeV, the 95% CL allowed ranges for the S and T 
parameters for U = are [19] 

Sbcst fit = 0.04 , a s = 0.09 , 

Tbest fit = 0.07, a T = 0.08, 

/W = 0.88, (3.1) 

where as,T denotes the la error and yO CO rr denotes the correlation between the two quan- 
tities. We discuss some regions of parameter space that are relevant for the remainder of 
this paper in the following. 

Obviously when all Yukawa couplings are zero the new particles decouple from the 
electroweak sector and AS = AT = 0, as long as the explicit mass terms are large enough 
for the S-T formalism to be a valid approximation. 
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The opposite limit, when only the Yukawa couplings are nonzero, is also well known. 
When the splitting between the members of the doublet goes to zero, i.e. = Y' c and 
Y£ = Y" , there is a positive contribution of to the S parameter of AS as 0.11, while 
AT = 0, which lies outside of the 95% CL region. Going away from the degenerate 
limit, AT grows proportional to the square of the mass difference of the members of each 
doublet, while AS receives corrections that are logarithmic in the mass difference. Due to 
the correlation between S and T this leads to a preference for a small splitting between 
the charged and neutral members of the doublets, with an upper limit on the splitting of 
about Am ~ 100 GeV if only one doublet is non-degenerate, as illustrated in the left plot 
of Fig 1. 

The degenerate case where the masses within one doublet are equal is also excluded, 
since here the positive shift in the S parameter is not compensated by a small positive shift 
in T. This region can be brought into agreement with the S parameter by increasing the 
Majorana masses, since they tend to reduce AS. 

Later we will be interested in scenarios with order one charged lepton Yukawa cou- 
plings, while keeping the neutral Yukawas small. In the absence of explicit vector-like 
masses such a scenario would be excluded by electroweak constraints. However here we 
can make use of the fact that the new sector decouples from the SM when the explicit mass 
terms rri£, m e and m v are increased while keeping the Yukawa couplings fixed. We therefore 
also expect that when m vec t or Y c v, the contributions to the S and T parameters should 
be reduced, since we are transitioning to the decoupling regime. This effect is shown in the 
right plot of Fig. 1 where in addition to the nonzero Yukawa couplings the explicit mass 
terms are set to mi = m e = m v = 400 GeV. For similar splittings between the charged and 
neutral Yukawa couplings, the contribution to the T parameter is suppressed compared to 
the chiral limit, such that a much larger region of parameter spaces is in agreement with 
precision constraints. 

In Fig. 2 the decoupling behavior is shown for S and T parameters separately. Here we 
choose a parameter region that will be relevant later, namely order one charged Yukawas 
with vanishing or small neutrino Yukawa couplings. To ensure that all neutrinos have 
sufficiently large masses the Majorana mass terms are set to M' = M" = 100 GeV. Further 
we choose Y' c = Y" = 0.8 and me = m e , so that the mass of the lightest charged lepton 
follows the simple relation 

M El = \m t - Y c v\ . (3.2) 

Only the region where m^ > Y c v is considered here. It can be seen in Fig. 2 that the S 
and T parameters decouple very fast in this regime. Sizable contributions only happen for 
charged lepton masses below 100 GeV, however even this region is still in agreement with 
the experimental limits due to the positive correlation between S and T. 

In addition to electroweak precision tests, there are also limits on the masses of new 
charged and neutral particles from the LEP experiments. The measurement of the width 
of the Z-boson at LEP constrains number of active neutrinos below 45 GeV to 3, which 
implies a lower bound of Mjy > 45 GeV for a Dirac and Mn > 39 GeV for an additional 
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Figure 2. Decoupling of S (red, solid) and T (blue, dashed) with increasing vector masses mg = m e , 
for Y' c = Y" = 0.8, Y^ = Y^ = and M' = M" = 100 GeV, as function of the lightest charged 
lepton mass m ei = mi — Y c v. Also shown is the ratio i? 77 (green, dotted) that is introduced in 
Sec. 4. 



Majorana neutrino [19]. Note however that these limits can be partially evaded when the 
coupling of the lightest neutrino to the Z boson is suppressed. 

Direct searches at LEP also impose lower bounds on masses of the new charged leptons. 
If we assume that the lightest new neutrino state N\ is stable on collider time-scales, the 
relevant decay is E\ — > W^*'N\ + where I denotes a SM lepton coming from the 
VF-boson decay. Production of EfE± pairs at LEP with subsequent decay to SM leptons 
and missing energy puts a limit of Me 1 > 101.9 GeV on the mass of the lightest new 
charged lepton, provided that AMi = Me 1 — MjVj > 15 GeV. When the mass difference 
between E\ and N\ is below this value, the resulting lepton can be very soft and therefore 
unobservable in the LEP detectors. In that case a weaker limit of Me 1 > 81.5 GeV applies 
for AMi > 8.4 GeV and M El > 63.5 for AMi > 7 GeV. No limits are available for 
smaller AMi, however limits from chargino searches [19] are available even for lower mass 
splittings, so that charged lepton masses below 80 GeV should be treated with care. 

In the absence of mixing between the new leptons and the SM leptons, there are 
no constraints from lepton flavor and lepton number violation, such that the electroweak 
precision tests and the direct search limits from LEP are the most important constraints 
on the parameter space for our model. 

4 Higgs Properties 

The couplings of the mass eigenstates to the Higgs boson are then given by the diagonal 
entries of the rotated Yukawa coupling matrix Ch = U^YcUn: 

C hll = Y^Ui u U Rll + Y?U* L21 U R21 , (4.1) 
C h2 2 = Y^U* L12 U R12 + Y?U* L22 U R22 . (4.2) 
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The new charged leptons, E% and E 2 contribute to the Higgs decay to photons at the one 
loop level. In the mass basis, the contributions to the decay are given by 



r 



OC 



where the loop functions for spin 1/2 and spin 1 particles are given by [20] 



Ai/ 2 (t) 
Mr) 



2(r + (r-l)/(r)) 
r 2 

o 3 3(2r-l)/(r) 



for spin - , 
for spin 1 , 



with 




itt + log 



i + yr 



l-VT 



for r < 1 , 
for r > 1 , 



(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



and m x is the mass of the particle running in the loop. It is instructive to consider 
the asymptotic values of the loop functions for Tj. < 1, i.e. when the new particles are 
much heavier than the Higgs boson. First, the SM contribution with = 125 GeV is 
Ai(t w ) = —8.3 from the W-boson loop and ^A 1 / 2 (Tt) = 1.8 from the top quark loop. 
Asymptotically A 1 / 2 (t -)• 0) = 4/3, while A 1 / 2 (r) > 4/3 for < r < 1. 

Since the new leptons don't affect the Higgs production channels, the effect on the 
di-photon search channel at the LHC is fully described by the ratio 



R 



77 



cr(pp — > h) T(h — > 77) 
<7sm(pp -> h) T(h ->■ 77) 



SM 



rfo ->• 77) 
r(/ l ^ 77 ) 



SM 



(4.8) 



In the limit of vanishing Dirac mass terms mi, m e , the pre factors Chuv/rrii in (4.3) go to 
one. It then follows that there is destructive interference between the dominant W boson 
contribution and the charged lepton loops. In this limit, we find a large suppression of the 
di-photon rate, i? 77 < 0.35, across the range of allowed Yukawa couplings. 

To better understand the conditions for obtaining an enhancement in i? 77 , we note that 
in the asymptotic limit the charged lepton contribution to the amplitude can be written as 



l l/2 



(o)£ 



Chuv 4 d 

= -v — logdet(A4) 

mi 3 dv 



4 2y c 'y c V 

3 y/y/'f 2 — mem e 



(4.9) 



where the first equality is a consequence of the Higgs low energy theorem [21-24]. Here 
it can also easily be understood by noting that = dA4/dv and remembering that 
tr(logAt) = log(det Ai). Constructive interference occurs when 4.9 becomes negative (and 
non-zero). For this to happen, both Yukawa couplings and both mass terms must be 
non-zero, and m^m e > Y/Y/'-u 2 must be satisfied. 
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Figure 3. Contours of constant i? 77 (green, solid) for Y' c = Y" = 0.8 as a function of the explicit 
mass terms mi and m e . The blue (grey) shaded region indicates a mass for the lightest charged 
lepton below 62.5 GeV, while the blue, dashed contours indicate a charged lepton mass of 62.5, 
100, 150, and 200 GeV. 

For order one Yukawa couplings, it is possible to obtain an enhancement of T(h — > 77) 
up to 50% compared to the SM prediction. This is illustrated in Fig. 3, where we choose 
Y' c = Y c " = 0.8 and vary the Dirac mass terms to determine the regions where the di-photon 
rate is enhanced or suppressed compared to the SM. 

Larger enhancements can be obtained in two different ways. First, one could increase 
the charged Yukawa couplings to values above one to get ratios i? 77 of two or larger. How- 
ever as we will see in the next section, such large Yukawas destabilize the Higgs potential 
below the TeV scale, such that an extension of the model would be required to realize such 
a scenario. Second, it would be possible to lower the mass of the lightest charged lepton 
below the LEP limit, and at the same time tune the mass of the lightest neutral state such 
that the decay E\ — > W^Ni produces a very soft lepton that escapes detection. In that 
case an enhancement of up to 70% can be obtained without further increasing the Yukawa 
couplings. 

To get a better idea about attainable values of i? 77 we have performed parameter 
scans with random choices for the mass terms and Yukawa couplings. In particular, we 
have varied rri£ and m e between GeV and 600 GeV and the Yukawa couplings Y' c and 
Y" in the range — 1.5. In Fig. 4 we show the correlation between the rate i2 77 and 
the mass of the lightest charged lepton. Obviously the maximal value of i? 77 increases 
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Figure 4. The ratio i? 77 for random values of m e G [0, 600] GeV, m e G [0, 600] GeV and F c ', F c " G 
[0,0.5] (blue), F c ', F c " G [0.5, 1.0] (red), and Y' c , F c " G [1.0, 1.5] (green). On the x-axis we show the 
mass of the lightest charged lepton. 

with the Yukawa couplings, however more importantly, it depends crucially on the mass 
of the lightest charged lepton. The di-photon rate increases for lower values of the lightest 
charged lepton mass, which is bounded from below by the LEP limits discussed in Sec. 3. 

Another interesting observation is that a small suppression of i? 77 can not be obtained 
in our model - either i? 77 > 1 or i? 77 < 0.3. The latter possibility, which includes the limit 
of rri£ = m e = 0, is excluded by the recent results from the LHC experiments at the 3<r 
level. Note however that it is possible to populate the region 0.3 < i? 77 < 1 when one of 
the explicit mass parameters is chosen to be negative. 

Experimentally the ratio i? 77 is not directly observable, since the strength of the 
photon signal also depends on the Higgs production cross section. Defining Rzz in analogy 
with (4.8), one should instead consider the ratio R^/ Rzz to establish a deviation of the 
Higgs branching ratios from the SM predictions. The prediction of our model is that the 
Higgs production cross section and the decays to vector bosons and SM fermions are not 
modified, and therefore R^/Rzz = -R77, in line with the current experimental data that 
shows no significant deviation of the signal rate in the ZZ* channel. 

A number of papers devoted to modifications of R^/Rzz through effects of new 
particles have appeared recently. General requirements on new physics to obtain a positive 
contribution are discussed in [24]. Analyses for specific models can be found in [25-30] for 
supersymmetric models 2 , in [31-37] for scalar extensions of the SM and in [38] for models 
with warped extra dimensions. 

2 In supersymmetric models it is also possible to obtain an enhanced i? 77 by suppressing the Higgs decay 
to bottom quarks [39-41]. Such scenarios lead to i? 77 /iizz ~ 1, which is a clear distinction from the 
predictions of our model. 
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Figure 5. One loop running of the gauge couplings (solid lines) compared with SM running (dashed 
lines). From top to bottom a^ 1 (blue), a^ 1 (red) and a^ 1 (green) are shown, and both green curves 
lie on top of each other. 

Since in all models the modifications are through loop effects, one possibility to distin- 
guish between different models is to compare i? 77 with Rz~/ since this Higgs decay is also 
loop induced [24] and eventually will be observable at the LHC [42]. Since the W-boson 
loop fully dominates the Higgs decay to Z*y, new fermions will in general have a smaller 
effect in this decay. Indeed we have calculated Rz-y in our model and find deviations from 
the SM prediction of at most 5% even in parameter regions where i? 77 is enhanced by 50% 
or more. 

5 UV stability and Unification 

The addition of particles charged under the electroweak gauge group and with nonzero 
Yukawa couplings will affect the running of gauge couplings and of the Higgs quartic 
coupling. At the one loop level, the RGE for the gauge couplings can be solved analytically 
and the running is given by 

ar\A) = ar\ mz ) - k log , (5.1) 

where a, = gf/(2ir) and 6i = f§, i>2 = -jp and 63 = —7, while the corresponding SM 
values are li, = ^ and —7. The evolution for scales up to A = 10 17 GeV is shown in Fig. 5. 
Compared with the standard model, the curves for \ are bent downwards, such that 
unification is slightly improved and happens at a lower scale, around A = 10 14 GeV. 

For the running of the Yukawa couplings, we neglect the light quarks and leptons, and 
also assume that the neutrino Yukawa couplings are negligible compared to the charged 
lepton Yukawa couplings. In this limit, the beta functions for the top Yukawa Yt and the 
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charged lepton Yukawas Y' c and Y" are given by [43, 44] 

Pt = ^Y {^-Y t 2 + (3Y t 2 + Y' c 2 + Yf) - 4vr f^ai + \a 2 + 8a 3 )) , (5.2) 

^ = 16^ + (3y ' 2 + ^ + ^"^ " An (f 01 + I° 2 )) ' (5 ' 3) 

where Yg_ = Y£, Y" . Finally the beta function for the Higgs quartic coupling, j3\ is 

0x = — ^ ( 12A2 " 4 ^ A ( 3 «i + 9 «2) + (4vr) 2 ^a 2 + ^ ai a 2 + ^a(\ + XG - h) , (5.4) 

G = 4 (3Y t 2 + Y c /2 + Y c " 2 ) , i? = 4 (3Y t 4 + Y c /4 + Y c " 4 ) . (5.5) 

Thresholds are treated naively by setting the top quark and lepton Yukawa couplings to 
zero below the respective particle masses. 

It is well known that for a light Higgs the quartic coupling at the weak scale is small, 
such that the RGE evolution is dominated by the large top Yukawa in the SM, and decreases 
as the RGE scale is increased. For = 125 GeV the Higgs quartic crosses zero around 
A = 10 9 GeV and remains at small negative values for higher RGE scales. While this poses 
a stability problem, are more refined analysis [45, 46] shows that the SM is stable against 
vacuum decay on time scales larger than the age of the universe. 

The additional Yukawa couplings in the present model make the model decrease faster. 
For the numerical study, we assume equal charged lepton Yukawa couplings and vanishing 
neutrino Yukawas, a parameter choice that is motivated by the results for the Higgs decay 
rates to photon pairs. The running of A for several choices of Y' c = = Y c is shown in 
Fig. 6. 

For Y c > 1 the running is fast enough to destabilize the Higgs potential at the TeV 
scale. This would imply the need for a UV completion at or below this scale, which could 
have sizable contributions to electroweak precision observables and Higgs branching ratios 
and thus question the validity of our analysis. To avoid this, we restrict our analysis to 
charged lepton Yukawa couplings smaller than unity. One should further note that small 
Yukawa couplings for the neutrinos, Y n < 0.5, will not have a significant effect on these 
results, since the destabilizing contributions to the beta function are quartic in the Yukawa 
couplings. 

The limits on the magnitude of the Yukawa couplings obtained here should be consid- 
ered conservative. An increased range of validity for the effective theory would be obtained 
if one allows metastability and demands instead that the vacuum decay time is longer than 
the age of the universe. 

6 Dark Matter 

So far, we have not considered the mixing of particles of the new sector with SM leptons. 
Such mixing can be prohibited by making all new leptons introduced in this model odd 
under a parity symmetry. A natural consequence then is that the lightest particle which is 
odd under this symmetry becomes stable on cosmological timescales, and, if it is produced 
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Figure 6. Evolution of the Higgs quartic coupling as a function of the scale A, for different values 
of the charged lepton Yukawa couplings (Yj, = Y" = Y c ), as indicated in the figure. Threshold 
were taken as 100 GeV, 173 GeV for the light charged lepton and top quark, respectively, while the 
threshold for the heavier charged lepton is given by the relation Me 2 = (2Y c v + 100 GeV), provided 
that mi = m e > Y c v. 

at some point in the early universe, will have a relic density that could be (part of) the 
dark matter of the universe. 

Leptonic dark matter candidates with unsuppressed couplings to the Z boson, such as 
ordinary fourth generation neutrinos, are excluded by limits from direct detection, since 
they predict a large interaction rate with nucleons [47]. Even when a nonzero Majorana 
mass is introduced to suppress the Z-boson mediated spin independent cross section, the 
limits on the spin dependent cross section still exclude Majorana neutrino dark matter 
with masses between 10 GeV and 2 TeV [47-50]. 

This constraint can be relaxed in the model considered here. For vanishing Yukawa 
couplings, the two singlet neutrinos or i/£ have no couplings to the Z boson. With 
nonzero Yukawas the singlet neutrinos will mix with the neutral doublet states and therefore 
can have nonzero but small couplings to the Z boson. For definiteness, we consider a 
scenario where only one neutral Yukawa coupling, Y^, is nonzero, while Y£ = m v = 0, so 
that IS sterile. The remaining neutrino states mix according to the 3x3 mass matrix 



with rows and columns corresponding to (r/£, uL, v'^). In the limit where the Yukawa 
coupling goes to zero, the two neutrinos from the SU(2) doublets form a Dirac neutrino M 
with mass me, while the singlet has a mass M' . 

As long as the Yukawa coupling is small, Y^v <C mg, and m! < mg, there will be two 




(6.1) 



heavy neutrino states N2 and N$ with masses Mjv : 



2.3 



mi and which are mostly composed 
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Figure 7. Left: Relic density of V\ as function of the mass m Vl — ttjdm and of the neutrino 
Yukawa coupling Y^. The region where the relic density is consistent with the observed value is 
shaded green (grey), the shape of the green band is explained in more detail in the text. The red 
lines indicate the old (solid) and new (dashed) limit on ergi from Xenon 100. Right: SI (black, solid) 
and SD (green, dashed) cross sections for scattering off protons in pico barn (pb). 



of v' L and z/o, while the third neutrino state iVi will have a mass ~ m! and will be 
dominantly singlet, with small mixings that generate a nonzero coupling to SM particles 
that are suppressed by powers of 



Y' 2 v 2 



1117 



m 



r2 



3.2) 



Motivated by the results of the previous sections, we fix the parameters of the charged 
lepton sector to mi = 205 GeV, m e = 300 GeV and Y' c = Y" = 0.8. In this case the 
lightest charged lepton will be close to the LEP limit, Me 1 = 105.8 GeV, which puts an 
upper limit on the mass of the dark matter candidate since -Mdm = Afzvi < Me 1 ■ 

A suppressed coupling to SM particles also implies that the annihilation rates are 
suppressed, such that the thermal relic density will be too large for very small Yukawa 
couplings. It is therefore important to accurately treat the regions where the annihilation 
is enhanced either resonantly or through co-annihilation effects. For this purpose we have 
implemented the model into the dark matter code MicrOmegas [51], and verified the results 
by comparing with analytic approximations [47, 52, 53] in regions where they are applicable. 
The couplings of the charged and neutral lepton mass eigenstates to the SM are given in 
the appendix. 

In Fig. 7 we show the relic density as a function of the dark matter mass todm and 
the Yukawa coupling Y' c . The two dips at Mdm ~ 45 GeV and Mom ~ 62 GeV are due 
to resonant annihilation into Z bosons or Higgs bosons respectively. Note that due to the 
narrowness of the Higgs boson the resonant suppression of the relic density only happens 
for Mdm slightly below m/j/2, where the resonant annihilation energy is obtained due to 
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the velocity distribution of the dark matter. Above Mdm = 80 GeV annihilation into 
VF+VF - pairs becomes kinematically accessible. Finally the strong suppression of the relic 
density for Mdm ~ 100 GeV is due to co-annihilation with the lightest charged lepton, 
which is not velocity suppressed. 

As mentioned above, for Majorana particles the spin independent (SI) scattering me- 
diated by Z-bosons is suppressed, such that the most important constraints come from 
the Higgs mediated SI cross section and from the spin dependent (SD) cross section 
from Z-boson exchange. For dark matter masses of order 100 GeV, the strongest limit 
on the SI cross section, crgi < (2 — 3) x 10 -45 cm 2 = (2 — 3) x 10 -9 pb, comes from 
the Xenon 100 experiment [54, 55], while the COUPP experiment [56] puts a limit of 
a SD,p < 5 x 10 -39 cm 2 = 5 x 10 -3 pb on the spin dependent cross section for scattering on 
protons. 

In Fig. 7 we show both the SI and SD cross sections for the same parameters as 
discussed before. While the limits on the SD cross section are not yet imposing a constraint 
on our model, the updated result of the Xenon 100 experiment [55] is sensitive to some of 
the models parameter space, thanks to an improvement over the previous result by about 
a factor of three. However one should keep in mind that astrophysical uncertainties and 
variations in the nucleon form factors can easily change these limits by up to a factor of 
two. 

A strong limit of <jsd < 5 x 10~ 40 cm 2 = 5 x 10~ 4 pb was reported by the CMS 
collaboration [57] from searches for mono-photons and missing energy. This limit is not 
directly applicable here since it assumes that the particles which mediate DM annihilation 
and scatterings have TeV scale masses. Interpreting these results in our model would 
require a detailed simulation of the signal in question, which goes beyond the scope of this 
paper. 

In regions where the annihilation rate is not enhanced, Yukawa couplings of order 
0.5 or larger are needed to obtain the correct relic density. This scenario is now strongly 
constrained by the most recent Xenon 100 results. In the light of the RGE evolution of 
the Higgs quartic coupling, a smaller Yukawa coupling is also much preferred. The region 
where Mdm < ^h/2, where the annihilation rate is enhanced resonantly, will allow the 
Higgs to decay to pairs of dark matter, leading to a sizable invisible Higgs width, which is 
disfavored by the most recent LHC results [7, 58] . 

The co-annihilation region seems to be most consistent with Higgs data and RGE 
constraints. Interestingly, this region also allows one to reduce the mass of the lightest 
charged lepton, since the strongest constraint, > 101.9 GeV only applies when the 

mass difference M El — is larger than 15 GeV. We have therefore analyzed this region 
in more detail, with special attention to the maximally attainable enhancement in the 
h — >■ 77 rate. We choose m£ = m e and Y' c = Y" = 0.8 since those choices maximize i? 77 . 
The mass of the lightest charged lepton is then simply given by 

M El = \mt - Y c v\ = \mt - 139.2 GeV| . (6.3) 
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Figure 8. LEP limits, relic density and i? 77 in the co-annihilation regime, as function of the 
lightest charged lepton mass m ei and Am = m ei — m Vl . The red shaded region is excluded by 
LEP, the green area is consistent with the relic density constraint for two choices of the neutral 
Yukawa coupling and the vertical lines are contours of constant i? 77 . See the text for details on the 
parameter choices. 

The results of our analysis are displayed in Fig. 8, in the Me 1 — AM plane, where 
AM = Me 1 —Mn-l . For a neutral Yukawa coupling of order 0.1 the preferred mass difference 
AM is between 5 GeV and 10 GeV, in which case most of the limits on additional charged 
leptons listed in [19] do not apply. An enhancement of the h — > 77 rate of more than 50% 
is now easily attainable in regions that are consistent with a thermal relic N\ . The regions 
outside of the green bands can be brought into agreement with the relic density constraint 
by modifying the neutral Yukawa coupling Y^. For a fixed coupling, the region below the 
green bands is allowed if other particles contribute to the dark matter relic density, since 
v\ is over-annihilated in these regions. 

7 Conclusions 

In this work, we have presented a compact extension of the SM by a vector-like family of 
leptons, and shown that it can provide a viable dark matter candidate and a source for 
the enhancement in the Higgs branching ratio to photon pairs, thus explaining two of the 
remaining mysteries of the SM. 




We have performed a detailed analysis of electroweak precision constraints, outlining 
the decoupling effects of the vector-like masses that suppress contributions to the S and 
T parameters even in the presence of sizable custodial symmetry breaking from Yukawa 
couplings. It is shown that the Higgs branching ratio can be enhanced if there is mixing 
in the charged lepton sector, while the case without mixing is excluded at the 3<r level by 
the observation of the Higgs boson in the di-photon channel. 

The running of the Higgs quartic coupling places an upper bound on the Yukawa 
couplings. If one demands that the model is stable up to energies accessible at the LHC, 
then Y c ', Y c " < 1, assuming that the effects of the neutral Yukawas are negligible. This leads 
to a strong correlation between the maximal attainable value for i? 77 and the mass of the 
lightest charged lepton, which has to be 0(100 GeV) in order to obtain a 50% enhancement. 
Such a light charged particle will eventually be in reach of the LHC experiments, thus 
making the model testable in the near future. 

In the absence of mixing between the new lepton sector and the SM leptons, the 
lightest neutrino is stable. To satisfy the relic density constraint while keeping the neutral 
Yukawa couplings small, either resonant annihilation or co-annihilation with the lightest 
charged lepton is needed. The latter possibility is preferred since it also avoids invisible 
Higgs decays to pairs of dark matter. The preferred mass difference AM = Me 1 — Mjvj 
for the co-annihilation scenario is between 5 GeV and 10 GeV. In this regime the strongest 
LEP limits on Me 1 do not apply, and charged lepton masses below 100 GeV and values of 
i? 77 ~ 2 become possible. 

At the LHC, the most promising signal will come from Drell-Yan production of E^E^ 
pairs that subsequently decay through a virtual W boson, E\ —¥ W*N\, giving rise to final 
states with leptons, jets and missing energy. In the co-annihilation regime the leptons and 
jets from the W* decay become increasingly soft, so that only events where sufficiently 
hard jets are radiated from the initial state partons will be detectable. In this case better 
sensitivity can be obtained from E1N2 and N1N2 production, followed by N2 —> WE\ or 
./V2 — > ZN\ decays where now the Z and W bosons can be on-shell. In addition the dark 
matter candidate can be searched for in direct detection experiments. The SI and SD 
cross sections are in reach of the next generation of experiments, but the search will be 
challenging in the co-annihilation regime. 

It remains to say that the model is mostly phenomenologically motivated, and should 
be viewed as an effective theory at the weak scale that will eventually be embedded in a 
more complete UV theory, which should also ensure vacuum stability beyond the 10 TeV 
scale. One option would certainly be a supersymmetric completion of the model with a 
somewhat higher supersymmetry breaking scale and superpartners above the TeV scale. 
Another intriguing possibility is the completion with a larger number of vector-like fami- 
lies of quarks and leptons that can give successful gauge coupling unification when three 
complete families are added to the SM [59]. 
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Appendix 

A Oblique Electroweak Parameters 

Lagrangian for masses of vector-like fourth generation neutrinos is 



Cum = \ (4 - R c < *Z) 



( Y£v m e \ 
Y^v m! m v 
m t Y/> 

\ m u Y£v m" j 



I 

R 
n 
R 

K7 



u 1 ' 



+ h.c. 



(A.l) 



Yukawa couplings and Majorana masses can be complex in general, making neutrino mass 
matrix a complex symmetric matrix. Such matrix can be diagonalized using Takagi de- 
composition as Mn D = V T M n V, where V is a unitary matrix. This gives mixing of left 
and right chiralities as follows 



/f l yA 




(<\ 


F L Y 2 


= V T 




F L Y 3 


< 


\¥ L N 4 ) 







and 



/PrA^A 






F R N 2 






IPrYs 











(A.2) 



Let the diagonal entries of M nD be Mn 1 _ 4 which are masses of Majorana neutrino states 
Ai_4 respectively. The Takagi decomposition produces a positive definite diagonal matrix 
which ensures the masses are positive. 

Lagrangian for masses of fourth generation charged leptons is 



■-eM 




+ h.c. 



(A.3) 



For chaged lepton states, mass matrix can be diagonalized by the singular value de- 
composition (SVD) given as M CD = U^ l M c Ur. The transformation between mass states 
Ei, E2 and flavor and chiral states e' L , e'l, e' R , after singular value decomposition is 
given as follows: 




and 



\WrE2 J 




(A.4) 
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The part of Lagrangian involving fourth generation couplings to and W$ boson is given 
as follows 

C n = f%fWfa + fZfWfa + fi%fWi<& + | 47^>r , (A.5) 

c n = g i^" w X - j^wfc + f <^wx - 9 i^ w >R ■ ( A -6) 

Substituting A. 2, A. 4 in the flavour basis Lagrangian given in A.5, we obtain the corre- 
sponding Lagrangian in mass basis as follows 

Similarly, we substitute A. 2, A. 4 in the flavour basis Lagrangian given in A. 6 to get 
£ W-3 =91 N rf 2 ^ * 

2 



52 £ %y 1J 2 2J fc w*s fc . (Ai 



Part of Lagrangian involving fourth generation couplings to photon is given as follows 
£ 7m = eAfiQ {Et-fEx + £27^2) = -eA„ {ErfEk + £37^2) . (A.9) 

From above described lagrangians, we can derive gauge boson self energies LTn (g 2 ) , II33 (g 2 ) 
and II3Q (g 2 ) as given in [16] and use the following formulae to find S and T. 



s = ^i [n 3 3( g 2 )-n 3Q (g 2 )] 



(A.10) 



q 2 =0 



T = (o) - n 33 (o)] . (A.11) 

C W b W IV1 Z 



We define 



A = M 2 2 x + M 1 2 (l-x)-x(l-x)< ? 2 , (A.12) 
b (Mi, M 2 , g 2 ) = jf log f dx , (A.13) 

b l (M 1 ,M 2 ,q 2 ) = [ xlog(£\ dx, (A. 14) 







b 2 (M 1 ,M 2 ,q 2 ) = J a: (1 - x) log dx = b 2 (M 2 , M u q 2 ) , (A.15) 

63 (M l5 M 2 , 0) = M|^(M 1 ,Af2,0) + M 2 6 1 (M 2 ,Af 1 ,0) 
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M*-M? + 2M 1 2 M|log (M 2 /A 2 ) - 2M 2 M 2 log (M|/A 2 ) 

/ 3 JMi, M 2 J = NlxM-i —3 , 

2 (M 2 - M 2 ) 

(A.17) 

where A 2 is an arbitrary regularization scale that will drop out of all physical observables. 
Note also that fz(M\,M\) = —1/6 is well defined in the limit M2 — > M\. Using these we 
can write 

2 

* S = E {\ U Lj\U Llk \ 2 + \U R2j \ 2 \U R2k \ 2 )b 2 (M Ej ,M Ek ,0) 
j,k=i 
2 

+ Y. Re { U ^U* Llk U R2j U R2k ) h (M Ej ,M Ek ) 
j,k=i 
4 

+ E (lni| 2 |^l fe | 2 + |^3i| 2 |^ fe | 2 )62(M^.,M A r fc) 0) 
j,k=i 
4 

j,k=i 
2 

- 2 E (l^i/ + I^T) h(M Ej ,M Ej ,0) + - , (A. 18) 

4tt4 c 2 M 2 T = 

2,4 

- 2 E (\U Llj \ 2 \V lk \ 2 + \U R2 .\ 2 \V 3k \ 2 )b 3 (M Nk ,M Ej ,0) 
j,k=i 
2,4 

+ 2 Re (f/ Ll3 £/£ 2 . FifeFsfc) M Ej M Nk b (M Ej , M Nk , 0) 
j,k=i 
4 

+ E (Vi;l 2 l^l 2 + \V 3] \ 2 \V 3 k\ 2 )b 3 {M Nj ,M Nk ,0) 
j,k=i 
4 

- E ^(Vi^ k V^ k )M N] M Nk b Q {M N]1 M Nk ^) 
j,k=i 

+ (I^Lnl 4 + \Ur 21 I 4 ) M 2 ^! (M £l ,M Sl ,0) + (|C/ Ll2 | 4 + |l^, 22 | 4 ) M E2 h (M E2 ,M E2 ,Q) 
+ {2\U Lll \ 2 \U R21 \ 2 + 2\U Ll2 \ 2 \U R22 \ 2 ) 63 (M El ,M E2 ,0) 
4 

- ^ Re (t^^t^yC^,) M Ej ME k bo [M Ej ,M Eh ,0) . (A.19) 
j,fc=i 

For the case where the explicit mass terms vanish, i.e. = U Rij = 5{j and Vj k 7^ only 
if j and k both belong to either {1,2} or {3,4}, A. 18 and A.19 reduce to the following 

vrS = - - 6 2 (M El , M El , 0) - 6a (M E2 , M Ei , 0) 
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2 4 

+ E l^'l 2 \ V ^\ 2 b2 {M Nj ,M Nk ,0) + l^f \V 3k \ 2 b 2 {M Nj ,M Nh ,0) , (A.20) 

j,k=l 3^=3 
2 2 

4^44 M 2 2 T = M% h (M Ej , M Ej , 0) - 2 £ | | 2 6 3 (M Wj , M Bl , 0) 
i=i i=l 

4 2 

-2 J] |F 3i | 2 & 3 (Mjv J ,M JSil) 0)+ ^ |y lj | 2 |y lfc | 2 6 3 (M 7 v J ,M JVfc ,0) 

3=3 j,k=l 
4 

+ E I y 3i 1 2 I ^ 1 2 ^3 (M^. , , 0) . (A.21) 

j,k=3 

A.20 and A.21 in turn reduce to the formulas given in [17] in the limit of vanishing 
Majorana masses. In the decoupling limit where all Yukawas and Majoranas go to zero, 
we get S = T = as expected for a vector-like generation as shown in Fig. 2. 
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